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Abstract 


Numerical viscous solutions based on an unstructured grid methodology are 
presented for a candidate high-speed civil transport configuration, designated as 
the Technology Concept Airplane (TCA), within the High-Speed Research (HSR) 
program. The numerical results are obtained on a representative TCA high-lift 
configuration that consisted of the fuselage and the wing, with deflected full-span 
leading-edge and trading-edge flaps. Typical on- and off-surface flow structures, 
computed at high-lift conditions appropriate for the takeoff and landing, indicated 
features that are generally plausible. Reasonable surface pressure correlations 
between the numerical results and the experimental data are obtained at free- 
stream Mach number = 0.25 and Reynolds number based on c R- = 8 x 10 6 
for moderate angles of attack of 9.7° and 13.5°. However, above and below this 
angle-of-attack range, the correlation between computed and measured pressure 
distributions starts to deteriorate over the examined angle-of-attack range. The 
predicted longitudinal aerodynamic characteristics are shown to correlate very 
well with existing experimental data across the examined angle-of-attack range. 
An excellent agreement is also obtained between the predicted lift-to-drag ratio 
and the experimental data over the examined range of flow conditions. 


Introduction 

In recent years there has been considerable 
progress made to advance the state-of-the-art unstruc- 
tured grid methodology for simulating nonlinear vis- 
cous flows over complete aircraft configurations. The 
primary interest for advancing such a technology is the 
inherent advantage and flexibility with which a tetra- 
hedral mesh can conform to a complex surface geome- 
try, along with the associated ease in generating the 
corresponding flow-field volume grids. In addition to 
the versatilities associated with the unstructured grid 
generation, complementary technology advancements 
in digital computers, both in tenns of speed and mem- 
ory, have also promoted notable progress in the 
unstructured grid algorithm developments for simulat- 
ing complex, three-dimensional, nonlinear viscous 
flows. Despite this encouraging progress, however, for 
various reasons, the application of this technology is 
not widespread among the general user community, in 
particular for viscous flow simulation. One may 
account for this lack of popularity because this is an 
emerging technology, and thus is not validated/ 
calibrated for complex flow characteristics often asso- 
ciated with realistic aerospace vehicles across flight 
speed regimes. As a result, the unstructured grid tech- 
nology is presently undergoing various stages of 
validation/calibration for different types of flows and 
configurations from fundamental (refs. 1 to 3) to com- 
plex (refs. 4 and 5) configuration problems. Such 


benchmark analyses are critical to building confidence 
and credibility in the numerical predictions for differ- 
ent classes of configurations and flow problems. 

One such method among various emerging 
unstructured grid technologies is the Tetrahedral 
Unstructured Software System (TetrUSS) (ref. 6) 
developed at Langley Research Center (LaRC). 
TetrUSS has four separate software packages: the flow 
solver USM3Dns, the grid generator VGRIDns, the 
surface geometry discretizer GRIDTOOL, and the 
solution postprocessor VPLOT3D. The Euler version 
of TetrUSS has been extensively validated and cali- 
brated for a variety of complex flows and configura- 
tions (refs. 7 to 11). In addition, a new adaptive 
unstructured grid Euler method has recently been 
developed (ref. 12) and coupled with TetrUSS. In ref- 
erence 12, the applicability of this new grid adaptive 
method has been successfully demonstrated to capture 
complex flow characteristics for a variety of aerospace 
configurations at flow conditions ranging from low 
speed to supersonic Mach number. In the past few 
years, however, the viscous capabilities of TetrUSS, in 
particular the flow solver USM3Dns and the grid gen- 
erator VGRIDns, have been maturing for simulating 
complex flows around geometrically complex config- 
urations (ref. 4). One main feature of the flow solver 
USM3Dns is the coupling of the Reynolds Average 
Navier-Stokes (RANS) equations with a wall function 
to simulate the boundary layer viscous flows. This 



approximation for viscous flow simulation is made 
primarily to alleviate the large computer memory 
requirements associated with solving the RANS equa- 
tions within the thin near-wall tetrahedral cells in the 
viscous sublayer region. The results from the applica- 
tion of this method for predicting the shock that 
induces separated flows at transonic high Reynolds 
number conditions for the ONERA M6 wing are pre- 
sented in reference 1 . In this report, the predicted sur- 
face pressure coefficients that are based on the wall 
function approximation are shown to correlate very 
well with experimental data, as well as with predic- 
tions from the full Navier-Stokes analysis. Further- 
more, the predicted surface flow patterns, computed 
with the wall function, are also reported to be compa- 
rable to those that are computed based on the full 
Navier-Stokes formulation. 

The study assesses the applicability of this method 
for predicting the complex separated flows associated 
with a High-Speed Civil Transport (HSCT) configura- 
tion, designated as the Technology Concept Airplane 
(TCA) within NASA’s High-Speed Research (HSR) 
program at low-speed, high-lift conditions. Such com- 
plex separated flows, emanating either from surface 
discontinuities (e.g., wing sharp leading edges and 
side edges of the deflected wing trailing-edge flaps) or 
from a smooth or round region of a surface (e.g., flap 
hinge lines and blunt wing leading edge) are consid- 
ered among the most challenging aerodynamic 
problems to simulate numerically. Although this 
particular method has not been validated/calibrated for 
such applications, attempts will be made in the 
solution analyses to demonstrate the method’s 

strengths and weaknesses in predicting some of the 
major flow physics associated with this class of 
aerodynamic problem. The accuracy assessments of 
the predicted results will be gauged through 
comparisons with the available experimental data 
obtained at LaRC. Figures in this report appear in 
color at http://techreports.larc.nasa.gov/ltrs/PDF/ 
1999/tp/NASA-99-tp209718.pdf. 

Symbols 

b reference span, 78.757 in. 

c reference chord, 56.972 in. 

C D drag coefficient, Drag/tq^ S ref) 


C L 

lift coefficient, Lift/iq^ S ref ) 

r 

'-'m 

pitching-moment coefficient, Pitching 
moment/Zq^ S ref c ) 

c p 

pressure coefficient, {p - p^/q^ 


free-stream Mach number 

P 

static pressure, lb/irr 

Poo 

free-stream static pressure, lb/in“ 

qoo 

free-stream dynamic pressure, lb/irr 

R c 

Reynolds number based on c 

^'ref 

reference area, 3060.0 hr 

* 

V 

wall-friction velocity, /t w / p , in/sec 

X,Y,Z 

geometry reference coordinate system 

y 

normal distance from viscous wall, in. 

y + 

inner-law variable, yv*/v 

a 

angle of attack, deg 

P 

viscosity, lb-sec/in 2 

v 

kinematic viscosity, (a/p, in 2 /sec 

P 

■2 

density, slug/in 

V 

wall shear stress, lb/irr 

Abbreviations: 

CFD 

computational fluid dynamics 

CFF 

Courant, Friedrichs, Fewy number 

CPU 

central processing unit 

FDS 

flux difference splitting 

FS 

fuselage station 

HSCT 

High-Speed Civil Transport 

HSR 

high-speed research 
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IGES 

Initial Graphics Exchange Specification 

LaRC 

Langley Research Center 

NAS 

Numerical Aerospace Simulation 

PVA 

primary vortex attachment 

RANS 

Reynolds Average Navier-Stokes 

SV 

secondary vortex 

SVSL 

secondary vortex separation line 

TCA 

Technology Concept Airplane 


TetrUSS Tetrahedral Unstructured Software System 

Experimental Data 

The wind tunnel experiment was performed in the 
Langley 14- by 22-Foot Subsonic Tunnel. The test 
was conducted with a 0.05-scale model of the TCA 
configuration, mounted on the floor of the tunnel, by 
using a stmt apparatus, as shown in figure 1. The pri- 
mary objective of the test was to obtain low-speed 
aerodynamic perfonnance data for the complete TCA 
configuration in the takeoff and landing flight modes, 
which are generally referred to as the high-lift 
conditions. In addition to testing the complete TCA 
model, a few experimental runs also were made with a 
more geometrically simplified version of the configu- 
ration, primarily to obtain data for computational fluid 
dynamics (CFD) method prediction assessments. This 
simplified geometry consisted of only the fuselage and 
the wing, incorporating different leading- and trading- 
edge flap deflections and excluding the engine 
components. 

In addition to the overall force and moment mea- 
surements, the model was instrumented for static pres- 
sure measurements over the wing upper and lower 
surfaces at various spanwise and chordwise stations on 
the starboard side. A schematic planform view of the 
TCA configuration, along with the spanwise pressure 
port locations to be used in the present analysis and 
correlations with computational predictions, are 
shown in figure 2. The measured longitudinal aerody- 
namic data were subsequently corrected to account for 
wind tunnel wall interference effects by using 
Fleyson’s wall correction method (ref. 13). Correc- 


tions were made also for the presence of the mounting- 
strut apparatus. The latter correction was based on the 
incremental effects of the measured forces and 
moments obtained in an earlier test conducted in the 
same tunnel on a 0.06-scale model of the same class 
vehicle, that is, the Boeing-designed Reference H 
configuration. 

The experimental data were analyzed prior to the 
development of the computational matrix to identify a 
candidate configuration and an appropriate range of 
angles of attack that are representative of the high-lift 
conditions for takeoff and landing. This analysis led to 
the selection of a simplified TCA configuration that 
consisted of the fuselage and wing, with fully 
deflected leading-edge flaps of 30° and fully deflected 
trailing-edge flaps of 10°. The numerical solid model 
representation of the selected configuration is shown 
in figure 3 from different perspective views. Subse- 
quently, the angles of attack of 5.7°, 9.7°, 13.5°, 17.2°, 
and 21° at = 0.25 and R- = 8 x 10 6 were selected 
for the computational analysis. 

Computational Approach and 
Attributes 

Flow Solver 

All numerical results presented in this report were 
computed by using the flow solver USM3Dns (ref. 1). 
The time-dependent, three-dimensional RANS equa- 
tions were solved to simulate the flow within the com- 
putational domain discretized by tetrahedral-mesh 
elements. USM3Dns is based on the cell centered, 
finite volume approach and uses an upwind-biased 
flux-difference-splitting (FDS) scheme (ref. 14) for 
spatial discretization of the primitive flow variables. 
Flow solutions are advanced by the implicit time inte- 
gration approach (ref. 15), with convergence acceler- 
ated to steady state by local time stepping. All 
computations presented in this report were obtained 
with the second-order FDS scheme, minmod limiter, 
and 10-stage subiterations for implicit Gauss-Seidel 
time stepping. 

For turbulent viscous flow simulation, the method 
solves the RANS equations along with the Spalart- 
Allmaras (ref. 16) one-equation turbulence model. The 
viscous flow solutions are then coupled with a 
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law-of-the-waU function, derived by Spalding (as pub- 
lished in ref. 17) to approximate the flow within the 
sublayer of the turbulent boundary layer (ref. 18). This 
boundary layer approximation provides a significant 
reduction in the overall memory requirement by elimi- 
nating the grids nonnally needed for resolving the 
boundary layer flow within the turbulent inner sub- 
layer (i.e., typically y + < ~20 to 50). Such an approxi- 
mation also improves the solution convergence 
characteristics caused by the reduced stretching of the 
near- the- wall cells, which generally adds stiffness to 
the numerical analysis. However, due to the inherent 
limitations of the wall function approximation, there 
are also genuine concerns about the accuracy of the 
numerical solutions in predicting the detailed flow 
physics, such as flow separation from a smooth sur- 
face. To predict this type of flow separation, which 
occurs because of the adverse pressure gradient within 
the boundary layer, requires the integration of the 
RANS equations to the configuration solid surface. As 
a result, efforts will be made in the present numerical 
analyses to examine the ability of the method in pre- 
dicting such a flow characteristic. 

Computational Grid Discretization and 
Attributes 

The wing full-span leading-edge flaps are 
deflected to 30°, and both trailing-edge flap segments 
(inboard and outboard) are deflected to 10° (fig. 3). 
This particular configuration, although geometrically 
simple, is thought to generate all major external flow 
physics associated with the takeoff and landing condi- 
tions experienced by the complete TCA vehicle. As a 
result, this baseline configuration, along with the cor- 
responding experimental data, is chosen for the 
present numerical analysis. The baseline TCA geome- 
try is defined parametrically for the wing and the fuse- 
lage in a fonnat known as Initial Graphics Exchange 
Specification (IGES), reference 19. This geometry 
definition served as the database for ah subsequent 
grid generation processes. GRIDTOOL (ref. 20) is pri- 
marily used to discretize the geometry into various 
surface patches which are then fed into VGRIDns 
(refs. 21 and 22) to generate the initial surface triangu- 
lation by the advancing front method. The initial sur- 
face triangles, also referred to as the initial advancing 
front, are then read back into the GRIDTOOL pro- 
gram for projection onto the initial database defined in 
IGES fonnat. The projection process of the initial 
front onto the initial database is required to preserve 


the proper surface curvature within the interior of a 
given patch. Finally, the projected surface triangles 
and the defined far-field boundary patches are used by 
VGRIDns to generate the complete volume grid (i.e., 
viscous and inviscid) within the computational 
domain, along with the corresponding grid and face 
connectivity files for the flow solver. The flow-field 
grid generated for the present analysis effectively used 
the newly developed feature of VGRIDns, allowing 
stretching of the grids in any direction on the surface 
or in the field. The present computational grid is 
stretched in the spanwise direction along the wing 
leading and trailing edges. The grid-stretching rate 
factor is strategically chosen for various regions of the 
configuration to efficiently model the local geometric 
features or the expected flow physics. The grid- 
stretching rate factor is varied anywhere between 1.05 
to 1.25 for different regions, and generally, the higher 
rate factor is applied to the field grid in the radial 
direction beyond the viscous advancing-front layers. 

The far-field boundary faces of the computational 
domain, which is shaped like a rectangular parallelepi- 
ped, are located at about 20 c upstream and 
downstream, 18c spanwise, and 9c above and below 
the numerical model (fig. 4(a)). For reference, 
c = 56.972 in., and the origin of the reference coordi- 
nate system (X = 0, Y = 0, and Z = 0) is defined to 
coincide with the fuselage nose apex. The overall con- 
figuration body length is about 195.6 in. in 0.05-scale 
model dimension. In the present numerical analysis, 
the external flow field is assumed to be symmetrical 
about the configuration plane of symmetry; thus, only 
half the configuration is modeled. Near-field and 
close-up views of the surface triangles and the grids in 
the plane of symmetry are shown in figures 4(b) 
through 4(f) from various vantage points. Suitable 
grids are distributed on the surface to capture either 
the various curvatures and discontinuities (e.g., flap 
hinge lines, trailing-edge-flap gaps, leading and trail- 
ing edges, and wing-fuselage juncture) and the 
expected flow characteristics. The computational 
model consisted of about 60000 surface boundary tri- 
angles, about 870000 tetrahedral cells within the vis- 
cous boundary layer region, and about 780000 cells in 
the inviscid computational domain. The computational 
grid contained about 18 to 21 tetrahedra cells within 
the boundary layer (see representative in the plane of 
symmetry, fig. 4(b)) with the first grid spacing sized 
nonnal to the surface to yield y + ~ 20 at approximately 
midfuselage body length. 


4 



Because of the wall function implementation, the 
USM3Dns slip velocity boundary condition is used on 
the configuration’s solid surfaces. Furthermore, the 
characteristic boundary condition is invoked on the 
computational domain far- field boundaries, along with 
the outflow plane. The free-stream boundary condition 
is specified for the inflow plane to initiate the 
computation. 

Numerical Results and Analyses 

The computations were initiated with the focus on 
the a = 13.5° case to develop an application procedure 
that could be applied to obtain consistent numerical 
results for the other cases. In the following first sec- 
tion, representative results will be discussed for solu- 
tion convergence, method perfonnance characteristics, 
and solution sensitivities (i.e., surface C p , integrated 
forces and moment) to the number of iterations. The 
predicted surface and off-surface flow characteristics 
will be presented in the second section for a typical 
solution. A more quantitative analysis and discussion 
on the predicted surface pressure coefficients, forces 
and moment, and corresponding correlations with 
experimental data will be included in the third and 
fourth sections. 

Solution Convergence and Sensitivity 
Characteristics 

All present computations were performed on the 
Numerical Aerospace Simulation (NAS) Cray-C90 
computer platform located at Ames Research Center. 
On this machine, the algorithm required about 44 psec 
per iteration per cell and about 240 million words of 
memory for a typical computation. A sample of solu- 
tion convergence characteristics is shown in figure 5 
for a = 13.5°. This figure shows the log of the total 
residuals (left vertical axis) and total integrated forces 
and moment (right vertical axis) as a function of itera- 
tions used to advance the solution. Figure 5 also shows 
the Courant, Friedrichs, Lewy (CFL) number strategy 
used to advance the solution and the total central pro- 
cessing unit (CPU) time requirement for the solution 
development. Computation for a typical solution is 
started from free-stream conditions with the CFL 
number initially set to 0. 1 and subsequently ramped up 
linearly in two stages to a final value of 150, as shown 
in figure 5. All computations are performed by using 


eight Cray processors. The required total CPU time for 
a typical computation is also shown in figure 5. 

The results shown in figure 5 indicate that a nomi- 
nal solution convergence apparently can be achieved 
with about 900 iterations, in which the total residuals 
are reduced about two orders of magnitude, and oscil- 
lations in the overall forces and moment are at negligi- 
ble levels. However, further analyses indicated that 
there are clear solution sensitivities on the computed 
surface pressure coefficient and on forces and moment 
with additional solution iterations. Figure 6 shows typ- 
ical solution sensitivities on computed surface pres- 
sure coefficients at eight fuselage stations (fig. 2) for 
three different iteration levels. Though small at this 
particular condition, the solution sensitivities become 
more pronounced at other angles of attack and are par- 
ticularly reflected in the computed longitudinal aero- 
dynamic characteristics shown in figure 7. The results 
shown in figure 7 clearly indicate the need for advanc- 
ing the solutions to at least 1800 iterations where the 
solution sensitivities to further iterations become neg- 
ligible. Therefore, a typical convergence characteristic 
plot, such as the one shown in figure 5, cannot be used 
solely to establish solution convergence; it is also nec- 
essary to examine at least one other figure of merit 
such as the surface pressure distribution, surface and 
off-surface flow characteristic, and surface skin fric- 
tion distribution. All subsequent results presented in 
this report will be based on final converged solutions 
with 2500 iterations. 

Typical Flow Features 

The surface pressure coefficient contours com- 
puted at a = 13.5°, M = 0.25, and R- = 8 X 10 6 are 
shown in figure 8 for the configuration upper and 
lower surfaces, as well as for the computational grid in 
the plane of symmetry. In addition, figure 8 shows the 
computed total pressure contours in various cross-flow 
planes that coincide with fuselage stations in which 
the surface pressure coefficients are measured (fig. 2), 
along with an additional plane just aft of the wing 
trailing edge (FS 3000). 

The flow over the wing upper surface generally 
can be characterized by the presence of two well- 
organized vortices (figs. 8(a) and (b)). The first vortex 
appears to form near the mid inboard wing leading 
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edge and is convected longitudinally downstream, 
roughly aft of the wing leading-edge crank. The sec- 
ond vortex forms over the outboard wing as the flow 
separates along the leading edge of the deflected flap. 
The low-pressure footprints associated with both 
vortices are clearly reflected on the surface pressure 
distribution contours. In addition to the two prominent 
wing vortices, there also appear to be two weaker vor- 
tices (i.e., lower total pressure viscous loss and higher 
surface pressure footprint) that form over the fuselage 
and roughly inboard of the inboard-wing-flap hinge 
line. The latter vortex originates from the wing apex 
region and appears to migrate downstream and coa- 
lesce with the stronger inboard wing vortex aft of the 
leading-edge crank near the wing trailing edge. The 
complex upper surface flow structures near the wing 
trailing edge will be discussed further in conjunction 
with the computed surface flow patterns. The flow 
over the configuration lower surface (figs. 8(c) and 
(d)) is characterized by a fairly benign C p distribution 
ranging from =0.0 to =0.2. In addition, the narrow 
band associated with the total cross-flow pressure con- 
tours over the wing lower surface can be attributed to 
the viscous losses within the boundary layer region. 

The simulated upper surface flow pattern, 
superimposed over the computed surface pressure 
contours, is presented in figure 9. Note the slight 
change, relative to previous figures, in the color map 
(shown in gray scale) used to display the surface 
pressure contour. See color version at http:// 
techreports.larc.nasa.gov/ltrs/PDF/1999/tp/NASA- 
99-tp209718.pdf. The surface flow traces are simu- 
lated by releasing particle seeds at the nodes of each 
surface triangle and restricting the tracing to the sur- 
face. The combination of surface C p and surface flow 
results are interesting from an analysis standpoint 
because the connection between the surface pressure 
footprints for the expansion-compression regions to 
the vortical-flow separation line (converging stream- 
lines) and attachment line (diverging streamlines) are 
highlighted. However, particular attention should be 
given to the interpretation of the predicted surface 
flow pattern caused by the inherent limitations associ- 
ated with the wall function approximation. Although 
the predicted surface flow pattern seems plausible, it is 
difficult to gauge the accuracy of the predicted surface 
flow pattern, mainly with respect to the flow separa- 
tion lines, due to the lack of complementary experi- 
mental data. With the present methodology, the 


prediction of a flow separation line becomes particu- 
larly intriguing when it occurs within the boundary 
layer caused by the adverse pressure gradient over a 
smooth surface. An example of such a viscous bound- 
ary layer flow separation is clearly demonstrated and 
discussed in the next paragraph. 

The same results as those presented in figure 9 for 
the surface C p and flow pattern are shown in the mid- 
dle part of figure 10 from an oblique rear vantage 
point, along with the computed total pressure contours 
in two cross-flow planes located near the wing trailing 
edge. No complementary experimental data exist for 
direct comparison with the predicted flow characteris- 
tics; however, the results shown in figure 10 are ana- 
lyzed solely to demonstrate the plausibility of the 
computed flow structures. The total pressure contours 
clearly show both the inboard and outboard wing pri- 
mary leading-edge vortices, along with the corre- 
sponding surface-pressure suction footprints (denoted 
by high-gradient regions) and the surface-flow stream- 
line pattern. The total pressure contours at both fuse- 
lage stations indicate the presence of a weak vortex off 
the side of the fuselage; however, the corresponding 
effects on the surface pressures and the surface flow 
pattern (i.e., the flow separation streamline) are not 
discernible. Similar flow analysis over the wing indi- 
cates flow characteristics that include the primary vor- 
tex attachment (PVA) lines for both the inboard and 
outboard wing (i.e., attachment lines are highlighted 
schematically with dotted arrow-lines for both the 
inboard and outboard wing vortices). Also, note the 
complex flow pattern around the wing tip region, a 
massive spanwise flow around the outboard wing trail- 
ing edge, and the formation of a small vortical flow 
over the gap between the wing and the inboard edge of 
the deflected outboard trailing-edge flap (see the total 
pressure contour plot at FS 2790). Furthermore, it is 
interesting to reveal that a secondary vortex (denoted 
by SV in the total pressure contour plots) appears to 
form under the inboard wing leading-edge primary 
vortex with a well-defined secondary vortex separa- 
tion line (denoted by SVSL in the figure). It is also 
important to notice that there is clearly no evidence of 
any local surface discontinuities that could trigger 
such a secondary vortex separation line. Although the 
prediction of such a flow characteristic with the 
present methodology (i.e., wall function coupled with 
the RANS equations) is intriguing; nevertheless, its 
manifestation is being reported here for the first time. 
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One plausible explanation for the presence of this sec- 
ondary vortex separation line can be attributed to the 
lack of a corresponding stagnation line due to the three 
dimensionality of the local flow field. As a result, the 
application of a slip boundary condition at the surface 
is sufficient to capture such a secondary vortex separa- 
tion line. 

The simulated lower surface flow pattern, super- 
imposed over the computed surface pressure contours, 
is presented in figure 11. The results generally indicate 
fairly benign and attached flow characteristics over the 
entire lower surface of the configuration. The lower 
surface pressure distribution also shows the expected 
flow compression near the deflected wing trailing- 
edge-flap hinge lines. 

Surface Pressure Prediction and Correlation 
With Data 

The computed surface pressure coefficients, along 
with the experimental measurements for a = 13.5°, 
M m = 0.25, and R- = 8 x 10 6 , are shown at eight fuse- 
lage stations in figure 12 (see fig. 2 for fuselage sta- 
tions). In an attempt to relate the pressure distribution 
to the configuration geometry, figure 12 also presents 
the corresponding geometrical cross sections for each 
respective fuselage station. Note that these geometric 
cross sections also graphically reveal the relative size 
of the inboard wing leading-edge radii (FS 1550 
through FS 2400), as compared with the sharp leading 
edge on the outboard wing (FS 2660 through 
FS 2790). In these pressure plots, the geometry dimen- 
sions for the local semispan and the Z-coordinate are 
given with respect to the 0.05-scale model. In addi- 
tion, no experimental surface pressure data were 
obtained on the fuselage; however, the numerical pre- 
dictions on the fuselage portion are presented for 
completeness. 

The computed upper surface pressure distributions 
at FS 1550 through FS 2010 are characterized by a 
narrow suction peak just aft of the wing leading edge, 
followed by a relatively mild flow expansion that 
appears to be at or just inboard of the wing-flap hinge 
line. The numerical results compare reasonably well 
with the experimental data at these three fuselage sta- 
tions. The computed upper surface pressure distribu- 
tions at FS 2240 and 2400 clearly indicate a wider 
band of low pressure that is associated with the 


inboard wing leading-edge vortex, which roughly cov- 
ers the entire flap local span. The leading-edge suction 
peak at these two fuselage stations is followed by a 
region of mild flow expansion, which appears to cor- 
respond to the wing apex vortex (fig. 8) inboard of the 
flap hinge line. The general character of the computed 
pressure distributions at FS 2240 and FS 2400 agrees 
with the experimental measurements; however, the 
magnitude of the suction peaks is underestimated for 
the wing primary leading-edge vortex. The computed 
surface pressure distributions over the last three fuse- 
lage stations also indicate a clear suction peak that is 
associated with the inboard wing primary leading- 
edge vortex, which appears to move slightly outboard 
at each consecutive fuselage station. The latter suction 
peaks are located at a local semispan of approximately 
21, 22, and 24 in. for FS 2510, 2660, and 2790, respec- 
tively. In general, the numerical predictions for these 
suction peaks do not correlate well with the experi- 
mental measurements. The computed upper surface 
pressure distributions that are associated with the out- 
board wing vortex over these last three fuselage sta- 
tions are predicted much better. 

The computed lower surface pressure distributions 
at all eight fuselage stations are predicted very well 
except at FS 2240 and 2400. At these two stations, the 
disagreements between the lower surface pressure pre- 
dictions and measurements are confined to the inboard 
section because of differences in the geometrical mod- 
eling. A possible source for this disagreement can be 
attributed to the geometrical differences between the 
numerical model and the wind tunnel model. The wind 
tunnel model consisted of a floor-mounted stmt appa- 
ratus (fig. 1) that was not included in the numerical 
model. Hence, the flow expansion around the stmt on 
the wing lower surface, causing lower surface pressure 
values, is not correctly simulated in the present com- 
putational study. 

Similar results for the computed surface pressure 
coefficients and the correlations with the experimental 
data are presented in figures 13 and 14 for a = 9.7° 
and a = 17.2°, respectively. These results are also 
based on the same Mach number and Reynolds num- 
ber shown in figure 12. At a = 9.7°, the surface pres- 
sure correlations between the predictions and the 
experimental data clearly show a much better agree- 
ment than the results presented previously for 
a =13.5° at all fuselage stations. However, at 
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a = 17.2° (fig. 14), the computed surface pressure dis- 
tributions deviate substantially from the experimental 
data at all fuselage stations. With increasing angle of 
attack, the complex surface flow separation from the 
round leading edges of the inboard wing and the 
subsequent formation of the primary vortex system are 
known to be very difficult to predict, even for a simple 
delta wing (for an example, see ref. 23). It should be 
noted, however, that the surface pressure disagree- 
ments do not appear to be as severe in the region 
associated with the outboard wing primary leading- 
edge vortex and the corresponding suction peaks (i.e., 
FS 2510, FS 2660, and FS 2790). The favorable pres- 
sure comparison in that region can perhaps be attrib- 
uted to the sharp leading-edge discontinuity, which 
fixes the flow separation line there, with the ensuing 
shear layers to fomi the outboard wing primary vortex. 

Force and Moment Prediction and Correlation 
With Data 

The computed longitudinal aerodynamic charac- 
teristics are presented in figure 15, along with the 
complementary experimental data obtained at 
M <*, = 0.25, and R- = 8 x 10 6 . In general, the correla- 
tions between the predictions and experimental data 
are very good across the range of angles of attack. The 
computed lift coefficients correlate very well with the 
measured data for moderate angles of attack; however, 
they appear to be overpredicted for a > 13.5°. It 
should be noted that the computed lift-oc curve exhibits 
a fairly linear slope, whereas the corresponding exper- 
imental data curve indicates a break in the slope near 
a = 5°, above which the slope is slightly less. The 
computed pitching moments and the drag ratios also 
indicate a very good agreement with the experimental 
data across the examined range of angles of attack, 
both in terms of magnitude and trends. 

The computed lift-to-drag coefficients are shown 
in figure 16, along with the corresponding experimen- 
tal data. The computational results clearly indicate an 
excellent agreement between the predictions and the 
measured data over the examined angle-of-attack 
range. This finding is particularly important because 
the CJC D ratio is often considered a figure-of-merit 
for assessing the high-lift aerodynamic performance of 
a given vehicle design. Within the HSR program, for 
example, the lift-to-drag ratio was used as a design 
parameter to be maximized, at a given C L (0.5 used for 


the TCA configuration), to arrive at appropriate con- 
trol surface deflections that provided the best aerody- 
namic performance for takeoff and landing conditions. 

Concluding Remarks 

Numerical viscous solutions based on an unstruc- 
tured grid methodology are presented for a candidate 
high-speed civil transport configuration, designated as 
the Technology Concept Airplane (TCA) within the 
High-Speed Research (HSR) program. The numerical 
results are based on the solutions to the Reynolds aver- 
age Navier-Stokes equations, coupled with a wall 
function, to approximate the flow within the sublayer 
of the turbulent boundary layer. The turbulent solu- 
tions, with the Spalart-Allmaras one-equation model, 
are obtained on a representative TCA high-lift config- 
uration that consists of the fuselage and the wing with 
deflected full-span leading- and trailing-edge flaps. In 
addition to the conventional approach for assessing 
solution convergence, the sensitivity of the computed 
surface pressures, with solution development, are also 
shown to be a good measure for establishing a 
solution-convergence level of completeness. 

Typical on- and off-surface How structures, com- 
puted on a representative high-lift configuration at 
conditions appropriate for the takeoff and landing, 
indicate features that are generally plausible. The abil- 
ity of the present numerical method to predict second- 
ary vortex separation-line-flow structure is reported 
for the first time. The prediction of this particular flow 
characteristic with a wall function approximation is 
found to be intriguing; hence, a possible explanation is 
presented for its formation. Because of the lack of 
complementary experimental data on the flow physics, 
the general plausibility of the predicted flow structures 
is primarily established by intuitive analysis and inter- 
pretation of the measured surface pressure distribu- 
tions. Such complementary experimental data on 
fundamental flow physics (e.g, surface flow pattern 
and off-surface velocity profiles) on this class of vehi- 
cles are essential ingredients used to gauge the accu- 
racy levels of the numerical predictions. 

Reasonable surface pressure correlations between 
the numerical results and the experimental data are 
obtained at = 0.25 and R- = 8 X 10 6 for moderate 
angles of attack of 9.7° and 13.5°. However, above 
and below these angles of attack, the correlation 


8 



between computed and measured pressure distribu- 
tions start to deteriorate over the range examined in 
the present study. The surface pressure disagreements 
are shown to occur in regions that were primarily asso- 
ciated with the vortical flows that emanate from the 
round leading edge of the inboard wing. However, the 
predicted longitudinal aerodynamic characteristics are 
shown to correlate very well with existing experimen- 
tal data across the examined angle-of- attack range. An 
excellent agreement is also obtained between the pre- 
dicted lift-to-drag ratio and the experimental data over 
the examined range of flow conditions. 

Finally, it is demonstrated that the present unstruc- 
tured grid methodology is a viable engineering tool for 
predicting the aerodynamic performance characteris- 
tics, with reasonable accuracy, for such vehicles at 
typical low-speed, high-lift conditions. However, it is 
recommended, as part of a future effort, to plan for a 
focused wind tunnel test to obtain a comprehensive set 
of on- and off-surface flow-field data, surface pressure 
measurements, and force and moment data, with suffi- 
cient detail and accuracy to validate/calibrate various 
computational codes. Such experimental data pres- 
ently do not exist for slender configurations such as 
the HSCT or the high-performance military class of 
vehicles at low speeds. 
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Figure 2. Schematic planform view of TCA and chordwise stations for experimental pressure measurements. 


Figure 1. The 0.05-scale TCA wind tunnel model in Langley 14- by 22-Foot Tunnel. 


Fuselage station (FS) 
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(b) Lower surface view. 

Figure 3. Solid model representation of selected TCA configuration. 


(a) Upper surface view. 
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(a) Overall computational far-field boundaries, 
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(b) Near-field view of grids on configuration upper surface and plane of symmetry, 
Figure 4. Computational grid for TCA high-lift configuration. 







(c) Oblique front view of upper surface grids on outboard wing. 



(d) Oblique rear view of upper surface grids on outboard wing. 
Figure 4. Continued. 
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(e) Near-field view of grids on configuration lower surface and plane of symmetry. 



(f) Oblique front view of grids on lower surface of outboard wing. 
Figure 4. Concluded. 
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Local semispan, in. Local semispan, in. 


(a) FS 1550 through FS 2240. 



Local semispan, in. Local semispan, in. 


(b) FS 2400 through FS 2790. 

Figure 6. Effects of solution advancement on computed surface C p \ 0C = 13.5°; = 0.25; R- = 8 X HZ’. 
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a, deg C D 


Figure 7. Effects of solution advancement on computed longitudinal characteristics. 
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(a) Upper surface view. 



(b) Aft-wing upper surface close-up view. 

Figure 8. Computed surface pressure and off-surface cross-flow total pressure contours; a = 13.5°; = 0.25; R- = 8 x 10 6 . 
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(c) Lower surface view. 


(d) Aft-wing lower surface close-up view. 
Figure 8. Concluded. 
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(a) Upper surface view. 



(b) Aft-wing upper surface close-up view. 

Figure 9. Computed surface pressure coefficients and flow pattern; a = 13.5°; M :Xj = 0.25; R- = 
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Figure 10. Computed on- and off-surface (total pressure) flow characteristics; OC = 13.5°; = 0.25; R- = 8 X H)' 1 . 
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(a) Overall view. 



(b) Outboard wing close-up view. 

Figure 11. Computed lower surface pressure coefficient and flow pattern; CX = 13.5°; = 0.25; R- 
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